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Abstract studies, most of which have been experimental, have 

been conducted to find the flow conditions that 


In the latest open-cycle gas-core concept, the 
large fuel-to-propellant density ratio of 10 or 
more, coupled with a vehicle acceleration of ap- 
proximately 0.02 "g" would cause a buoyancy force 
which would tend to reduce the amount of fuel con- 
tained in the cavity. This effect of buoyancy on 
fuel containment would depend on various engine 
flow parameters. It is not clear how these param- 
eters interact. Experimental flow tests covering a 
sufficiently wide range of parameters would be too 
costly. 

The piirpose of this analysis is to determine 
scaling laws for the buoyancy effect so that ex- 
perimental conditions can be related to the engine 
conditions. The fuel volxme fraction (proportional 
to the contained fuel mass) in a short coaxial flow 
cavity is calculated with a programmed numerical 
solution of the steady Navier-Stokes equations for 
isothermal, variable density fluid mixing. The so- 
lution includes buoyancy, pressure gradient, and 
viscous shear terms, as well as the effects of side 
and end walls. Solutions are obtained at a 
propellant-to-fuel flow ratio of 50 and a Reynolds 
number of 1000, for fuel-to-propellant density ra- 
tios ranging from 1.0 to 10, and for variofis ve- 
hicle accelerations. The main limitation of the 
analysis is that the numerical solution fails to 
converge for Reynolds numbers above about 1000. 

A dimensionless parameter B called the 
Buoyancy number was found to correlate the fuel 
volume fraction for large accelerations and vari- 
ous density ratios. This parameter has the value 
B = 0 for zero acceleration, and B = 350 for 
typical engine conditions. Calculation results 
show that for B less than 5.0 the fuel volume 
fraction, which depends strongly on density ratio, 
is independent of acceleration. However, for B 
greater than about -40, the fuel volume fraction is 
nearly independent of density ratio and decreases 
with increasing Buoyancy number. At a particular 
engine value of B = 350 the contained fuel mass 
is decreased by a factor of about 4.0 due to ve- 
hicle acceleration. 

Introduction 

The gas -core nuclear rocket is a proposed 
space propulsion system capable of high specific 
Impulse (1500 to 5000 sec) and a wide range of 
thrust. In the latest open cycle concept, de- 
scribed by Ragsdale^^-', shown in Fig. 1, the heavy 
uranim .fuel vapor fills the central region of an 
approximately spherical cavity and is STorrounded 
by lighter, faster moving hydrogen propellant. The 
propellant stream is heated to about 20,000° K by 
thermal radiation from the fissioning fuel core 
which is at a temperat\ire of about 60,000° K. 

One characteristic of any open-cycle concept 
is that some unfissioned fuel is lost. This occurs 
because the two gases are in direct contact inside 
the reactor cavity. Numerous fluid mechanics ; 


minimize fuel loss rates and maximize , the. contained ; 
fuel mass ( or equivalently the fuel vql'ume;::;,.f rac- 
tion). Most, recently the experiniental,),worKv:6£::,^;- 
Bennet and Johnson(2) has shown tha|:.',fuei:,'VoI'ito 
fractions of 50 percent at propellanttto-fue.I.i.tahS 
flow ratios above 50 can be obtained:.ih .coaXiSl.T.flow 
cavities. 

The recent emphasis on relatively. low,, thrust, , 
very high specific impulse, engines for Mars mis - 
sions(lJ calls for engine flow cqhditions with: very 
low average fuel velocities (about 0.01 ft/sec), and 
a fuel density that is about ten to twenty .times the 
average propellant density. Vehicle acceleration 
would be approximately 0.02 times earth gravity. 
These flow conditions can lead to large buoyancy 
forces on the fuel with the result that fuel con- 
tainment may be reduced by fuel acceleration out of 
the cavity. Such a "buoyancy effect" observed 
experimentally by Bennet and Johnson^''"-' in a two- 
fluid coaxial flow cavity. Although E -was rot a 
strong effect in the experiment, these invesiig-v.ors 
concluded that at conditions close to 'be ei the 
effect could be much stronger because ’ he rie-1 . ' - 
locities in the engine are much, less ' nan In Uit 
experiment . 

The effect of buoyancy on fuel containment 
would depend on many engine parameters, such as in- 
let velocity and density, cavity geometry, Raynctda 
number, and vehicle acceleration. It is not clear 
how these parameters interact. A thorough experi- 
mental study covering the full range of parametei-s 
would be much too costly, if not impossible, 'El.ore- 
fore analytical predictions are needed to generalize 
the experimental results so that one can extrapoi-uto 
the present observations to the engine conditions. 
This report describes an analysis that determines 
the functional dependence of contained fuel mass (or 
fuel volume fraction) on the fuel-to-propellant den- 
sity ratio and vehicle acceleration for the simpli- 
fied case of isothermal flow in a short cylindrical 
cavity. 

The purpose of this report is to analytically : 
determine the proper scaling laws for the "buoyancy 
effect" so that experimental conditions can be. .re- 
lated to engine conditions. Previous analyses :snch 
as for the ducted jet by Ghia(^) the free /jet:..^ 
Putrei'^) or the steady Navier-Stpkes equation ..soLu- : 
tion of the ducted jet by Shavit(?). were ruled .out 
for this study, either because the .equation: sets 
used neglected pressure gradients and recirculati'ng 
flows, or because the influence of the end walls 
(i.e., the nozzle constriction) was not included. 

It was felt that these factors might be important 
elements of the buoyancy effect. The .main: assump- : . 
tion of the present analysis is that the flow is 
steady, isothermal, and laminar. The analysis i’n- 
cludes a complete fluid force balance, including 
pressure terms, viscous shear terms, buoyanc5r 
terms, and effects due to a downstream end wall. 

The nmber of parametric calculaLlons :L.- ni. ■- 
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eGsarily ].imited because of the relatively long 
computing times for the ntimerical solution (typi- 
cally 15 minutes per case), which was due mainly to. 
numerical convergence requirements. Therefore, the.; 
following quantities were fixed as shown in Fig. 2; 
cavity L/D = 1, propellant-to-fuel mass flow ratio 
of 50, Reynolds number of 1000, and nozzle area 
ratio of 4 to 1. The generalized inlet velocity 
and fuel concentration profiles are also fixed. 

The fuel~to -propellant density ratio is varied from 
1.0 to 10.0 (and similarly the velocity ratio). 

As will be described in subsequent sections of this 
report, the buoyancy force can be characterized by 
a dimensionless Buoyancy number, B, which is pro- 
portional to the square root of vehicle accelera- 
tion. The B value in the calculations will be 
varied from zero to the largest value for which 
the program converges, or about B = 100. The com- 
puted results wi.ll be extrapolated to B values 
near 350, which correspond to the "0.02 g" level 
of the engine. As is shown in Fig. 2, except for 
Reynolds number, these flow conditions compare 
closely to the experimental conditions of Ref. 2, 
so that the results of this study can be applied 
to the coaxial flow experiment. 
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Symbol List 

generalized coefficients in the three 
flow equations 

generalized coefficients in the substitu- 
tion equations 

Buoyancy number, defined in equation (15) 
cavity diameter 

■ generalized dependent flow variable in 
the three flow equations 

smoothing function for inlet velocity and 
mass fraction profile 

vehicle acceleration 

length of cavity 

mass flow rate 

normal distance from wall 

cavity radius 


Reynolds number 


2(mp + iDp) 




inner and outer velocity buffer radius 

inner and outer mass fraction buffer 
radius 

radial co-ordinate 
Schmidt number 

specified axial velocity at the inlet 
axial velocity 

fuel volume fraction, defined as the 
fraction of the cavity filled with 
pure fuel if it were gathered together 
into a central volume 


y radial velocity 

X , axial coordinate 

y fuel mass fraction, local fuel density/ 

local mixture density 

p density 

(j. viscosity 

Q : : vorticity 

ijf stream fionction 

cp generalized dependent flow variable in the 

substitution equations 

Subscripts: 

F fuel 

i index for central node 

d,e,w,n,s indices for surrounding nodes 

K index for each of three flow equations 

P propellant 

0,1 wall node, adjacent inside node 

Analysis 

The cavity flow model to be analyzed is shown 
in Fig. 3. The analysis is restricted to a rec- 
tangular cylinder geometry as shown. The assump- 
tions of the analysis are: 

1. Steady, viscous flow, constant viscosity 

2. Isothermal flow 

3. Inlet and outlet flows are purely axial 

4. Impermeable, no-slip side and end walls 

Equation Set 

The system of equations to be solved is a re- 
written form of the variable density Navier-Stokes 
equations, the continuity equation, and the mass 
diffusion equation. The solution method is basic- 
ally the same one described in detail by Gosman, 
et al.v®). For sake of brevity only the main points 
of the analysis are repeated here. The rewritten 
equation set takes the form of three simultaneous, 
and similar, second order nonlinear partial differ- 
ential equations in the variables vorticity, stream 
function, and fuel mass fraction (concentration). 
These equations are generalized into the form 
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The first line in this general equation has the 
convective terms, the second line has the diffusive 
terms, and the third line has the general source 
tern including the acceleration body force term. 

The coefficients ajj, bjj, cjj, and dg- can depend on 
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x,r and the primary variable Fj^, and axe given in 
table I. The Schmidt number is taken as Sc = 1.51. 

The value of axial acceleration g in the body 
force term of d3_ (see table l) is specified in 
terns of the dimensionless Buoyancy number 
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(PF - Pp)gD 
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where the form of B is discussed later. 

Additional equations are for vorticity and 
stream function in terms of velocities, and for 
density in terms of mass fraction. 


u = ^ 

(2) 

pr Sr 

V = -A|t 

pr ox 

(3) 

/ Sv Su\ 

(4) 

PpPp 

~ pp + y(pp - pp) 

^ (5) 


The solution of this equation set requires 


that the values of Fj^ or SFjy'Sn 
all boundaries of the flow region, 
boundary conditions are: 


be specified on 
The centerline 


Sr I 


0, ij, = 0, = 0; at r = 0 ; (6) 

The outer no-flow wall boundary conditions are: 


= const, 
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where the wall and inside node are given the sub- 
scripts 0 and 1, respectively. The outlet 
boundary conditions are, for the outlet far from 
the corner: 


^3 


0 , = 0 , . 

ax Sx 


( 8 ) 


The boundary conditions of Eqs. (6) to (8) are 
quite general in that the same equations apply for 
any flow in the cavity. On the other hand the flow 
is strongly Influenced by the inlet bovmdary con- 
ditions. Therefore, these were chosen to resemble 
the inlet conditions in the experiment of Ref. 2. 
The iiilet velocity profiles and mass fraction pro- 
files were specified as smooth profiles with buffer 
regions (see fig. 3) and were numerically inte- 
grated to get the inlet stream function values. 
These are specified in the forms 

r < RYl, y = l.Oj RYl < r < EY2, y = 1 - f(r), 

r > RY2, y = 0 (9) 

r < EUl, u = Uj,; RUl < r < RU2, 

u = Uj, + (Up - Uj,)f(r), r > RU2, u = Up (10) 


il/(x = 0) 




pur cr 


( 11 ) 


where f(r) is a double parabola siiiOobLiiug I’lmitiou 
that runs from 0 to 1. Based on preliiiiinary C'l- 
culations, the buffer radii val_ues tcere chosen as 
EYl = 7/21, RY2 = 11/21, EUl = 10/21, 


RU2 X8/'21. 


Realistic vorticity inlet conditions cthat are 
compatible with, the above inJ.et conditions must 
also be specified. The inlet vorticity ..is. related 
to the "straightness" of the inlet :..flow through .t v 
Eqs. (2) to (4). The inlet vorticity was calcula- . 
ted from the expanded equation c 
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r r 
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where the term with S\];/5x is identically set^ to ■ 
0 to represent a purely coaxial , flowb at . :.the..:inletv . 
As will be discussed later, the .inlet iCondil'id^ 
are completely specified when RUl, RU2, RYl;,:(HYl 
Pp, Up, Pp, and Up are specified. : 

The fuel volume fraction is defined to be the 
fraction of the cavity filled mth pure fuel if it 
were gathered together into a central volume. 

This quantity is calculated from 
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numerical Solution Method 

The numerical solution is baseo on the mei 
detailed in Eef. 6. In simple terms, the .this. 
b,c) are written in terms of finite hiffereiu'.e, 
node points of a rectangular mesh covering the 
flow region. Central differences s,ro used for 
terms except the convective terms whtch use up' 
differences. (Upwind differenciiig improves th. 
stability of the numerical solution. ) The ros 
ing equation set can be written in the forju 
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where 
or y. 


q^i 

And 


■cpi - I '“i.tdl ‘ "1 

a=e,w,n,s 

is_one of the primary variables .ot/r, 1]/ 




or bj_ may depend on : q)r 


The upwind differences make the coefficient 
matrix of this equation se;t diagonally ^ dominant . :. . 
Thus if the coefficients a. i, and were con- 
stants (linear equation) ,. the solution by. vsucces'- . 
sive substitution (also IlmoWi as Gauss Seidel 
iteration) would be absolutely convergent. These 
coefficients will depend on local gradients, es- 


pecially velocity and density gradients 


Hence, 


.convergence may be difficult for those cases .with:., 
very large gradients. 

The stability of the equations (tendency to 
converge) were further improved by : several: methods. 
First, the buffer regions in the inlet veloc:ity : v 
and fuel concentration profiles were made as large 
as possible to reduce the gradients at the inlet. 
Secondly, the mesh was stretched in ine x direc- 
tion with node points closer at t.lie inlet than at 
the outlet. Thirdly, imderrelaxa bxon vuc nrv''fl i i 
dampen fluctuations in the cccf Cl ch-ntr, bM.Vi 'c . 
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iterations. After extensive trial runs it was 
i’oi,md that, for convergence over the large range of 
density ratios and velocity ratios to be studied, 
the stream function, local density, and vorticity 
source term dj^, req.uired tuaderrelaxation. ' 

Computer fTograni 

The numerical solution was programmed for. use 
on the IBM 7090/94 computer. The program, was 
based on the basic code described in Gosman, 
et al. (6). Substantial modifications and addi- 
tions ive;re made to this basic code. These in- 
cluded a subroutine for specifying and integrat- 
ing the inlet velocity, mass fraction, and stream 
function profilesj a subroutine for calculating 
the som'ce term including the gravity term 

for the vorticity equation,- a boundary condition 
subrouti,;;iej a velocity and density distribution 
subroutine I as well as input -output routines such 
as for writing tapes to restart calculations and 
to generate microfilm contour plots of stream 
fune.tion and fuel mass fraction. 

Calculations were usually initialized by 
starting: from: a previous solution on tape. This 
method reduced computing time by factors of 2 or 
3.: The iterations were run until the change in 

: between: successive, iterations was less than 
0.1 percent, A typical case calculation with a 
22 X 2:5 mesh required about 300 iterations with 
about 15 minutes of computing time. 

Results and Discussion 

The results of preliminary runs for deter- 
..mining.:. .the:., range of input parameters for oonverg-: 
enee, :ar:d the selection of inlet conditions are 
discussed : first. Then the results of parametric 
runs : at various density ratios and vehicle accel- 
erations: are discussed and analyzed. 

Limits: on: Range of Input Parameters 
for Convergence 

As: was discussed in the Analysis section, the 
numerical solution fails to converge when the 
equations: become strongly nonlinear because of 
large gradients in the flow variables. The condi- 
tions... for improved stability given in the Analysis 
section^ were used to obtain convergence for a 
fairly large range of parameters. Trial runs 
showed that the smallest flow ratios gave converg- 
ence at the highest Reynolds numbers. Thus, for 
this report the flow ratio was fixed at the 
smallest realistic value of propellant -to-fuel- 
flow ratio of 50. This gave convergence to 
Reynolds number as large as 10,000 with a density 
i-atio Ox unity. For the parametric calculations 
using fuel-to-propellant density ratios from 1.0 
to 10.0, a constant Reynolds number of 1000 was 
used to assure convergence. The runs for various 
buoyancy forces failed to converge for Buoyancy 
numbers, B, greater than 100. 

Effect of Inlet Conditions 

.Preliminary runs showed that for a fixed flow 
ratio and density ratio the cavity flow is 
strongly dependent on the inlet velocity and con- 
centration profiles. The object of these prelim- 
.inary runs was to determine those buffer radii 
for Eqs. (8) to (lO) which resulted in the largest 


fuel volimie fraction with good convergence. It 
was found that step profiles resulted in strongest 
recirculation and low fuel volume fractions. Pro- 
files with wide buffer layers resulted in least 
recirculation and increased fuel volume fractions. 
(In general strong recirculation coincided with 
low fuel volume fractions . ) The buffer radii val- 
ues which resulted in the largest fuel volume frac- 
tion values for fuel-to-propellant density ratios 
of 1.0 to 2.0, are given in the Analysis section. 
(These buffer radii were not re-optlmized for other 
density ratios, or for the nonzero buoyancy force 
cases.) It is worth noting that these analytically 
determined buffer values give about the same opti- 
mm inlet profiles as were experimentally deter- 
mined by Johnson^’^^. 

Streamline and Fuel Concentration Elots 

The effect of various density ratios and ve- 
hicle accelerations on the general cavity flow is 
best illustrated by the streamline and fuel con- 
centration contour plots in Figs. 4 to 6. These 
figures are from computer plots of the converged 
numerical solutions. The streamline plots in 
Fig. 4 and the fuel mass fraction contour plots in 
Fig. 5 show the large recirculation cell that is 
formed when the density ratio is increased from 
Pj,/Pp = 1.0 to 4.0 with no buoyancy force. This 
recirculation cell, clearly outlined by the zero 
streaniline in Fig. 4(b), convects a substantial 
amount of fuel away from the fuel plume, thus re- 
ducing the fuel mass in the cavity. This explains 
the reverse cxnrvature of the 50 percent contour in 
Fig. 5(b) as compared to Fig. 5(a). The fuel vol- 
ume fraction decreases from 0.072 to 0.034 as the 
density ratio increases from 1.0 to 4.0. 

The concentration contour plots in Fig. 6 
show the effect of vehicle acceleration on the 
shape and size of the fuel plume for a density 
ratio of Pp/Pp =2.0. Figure 6 shows that, as 
vehicle acceleration Increases from B = 0 to 65, 
the fuel plume is stretched out and gets much nar- 
rower. The fuel volume fraction decreases from 
0.059 to 0.0126, by nearly a factor of five. This 
fuel region stretching is the result of buoyancy 
force accelerating the fuel toward the downstream 
end. The concentration contours at B values 
near 100 for density ratios greater than 1.0 have 
shapes similar to those in Fig. 6(b). These flows 
are characterized by an absence of the recircula- 
ting cell mentioned above (with zero acceleration), 
and by substantially lower fuel volume fractions. 

It shoiILd be noted that the fuel volime frac- 
tions calculated for all cases of this report are 
factors of nearly 5 below the measured values from 
the coaxial flow experiments of Ref. 2. As dis- 
cussed above, the specified inlet boundary condi- 
tions were optimized by preliminary runs, to give 
the largest calculated fuel voi™e fractions. 

Thus no significant increase in calculated fuel 
voliane fractions is expected for different inlet 
conditions. Therefore it appears that the low cal- 
oiiLated values of fuel volume fractions could pos- 
sibly be the result of the calculations being 
limited to Reynolds numbers near 1000. This is 
somewhat confirmed by the fact that additional cal- 
culations for a density ratio of 1.0 at Reynolds 
number of 10,000 resulted in a fuel volume frac- 
tion of 0.108 compared to the value of 0.072 at 
the Reynolds number of 1000. The fact that the 
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analysis does not converge for very high Reynolds 
numbers, and that it results in consistently low 
fuel volume fractions, are the major shortcomings 
of this analysis. However, this investigation- is 
concerned mainly with the functional dependence of 
fuel volume fraction on various parameters. Thus 
for purposes of this report we are not as concerned 
with the precise values of fuel volume fraction as 
with the trends in the values. 

The variation of the flow with- mass: flow ra- 
tio and outlet diameter is also noted here. Based 
on a few calculations, the larger mass flow ratios, 
mp/mp, require larger values of velocity ratio, 
Up/Up, which results in stronger recirculation and 
correspondingly lower fuel volume fraction. The 
outlet diameter for the parametric runs was 0.48 
times the cavity diameter. When the outlet diam- 
eter was decreased to as small as 0.14 times the 
cavity diameter, the fuel volume fraction for a 
fixed engine acceleration was not significantly 
changed. 

Effect of Density Ratio on Fuel Volume 
Fraction at "Zero-g" 

The flow equations wre solved for fuel-to- 
propellant density ratios of 1.0, 2.0, 4.0, 
and 10.0. The flow ratio was fixed at mp/mp = 50, 
and the Reynolds number was 1000. The buoyancy 
force for these cases was assumed to be zero. The 
propellant to fuel flow ratio was adjusted over the 
range Up/Up = 16.3 to 183 to maintain the fixed 
flow ratio. The variation of fuel volume fraction 
with density ratio is shown in Fig. 7. The fuel 
volume fraction decreases from 0.072 to 0.017 as 
density ratio increases from Pp/pp = 1.0 to 10.0. 
The fuel volume fraction first decreases slowly as 
density ratio increases from 1.0 to 2.0, and de- 
creases faster (as the -3/4 power of density ratio) 
for density ratios beyond about 2.0. The two- 
region behavior is probably due to the fact that 
for the flow ratio of 50 recirculation starts at 
about Pp/Pp = :2 -0 and becomes stronger with in- 
creasing density ratio. The -3/4 power dropoff 
with density ratio is stronger than was calculated 
in Ref. 4 with a boundary layer solution (which 
did not include recirculation) where the variation 
was as the -1/2 power. The experimental results of 
Ref. 2 also showed about a -1/2 power variation, 
which was probably due to the fact that at the high 
Reynolds nurabers recirculation was weaker than in 
the present calculations. 



This parameter is called the Biojiiic, ■ iiumbo-- Li 
this report. The combination ins,dc ’J, u i .. ilc. L 
is a measure of the buoyancy Foj-ee tr ulic intr u 
force on the fuel. The square loot is aset l‘oi 
ease In plotting the B values wliicb ]’xnoo Or oi 
0 to about 350 . 

The fuel volume fraction is ploeLcd in 3'‘i;_ , 6 
as a function of Buoyancy number for various dc u- 
sity ratios. , The calculations range from B = i 
to about 100, where the upper ^^alue is the "urnicti- 
cal convergence limit. For a fixed density ratvio 
the fuel volume fraction is constant for small B 
values, and decreases as a constant power of B 
for large B values. Of particular importance is 
the fact that, for the combination of variables in 
B, the curves for the various density ratios: fall 
close together at large B values, which confirms 
the form of Eq. (15), 

At Buoyancy numbers above about 40, the fuel 
volume fractions are close enough to a single ex- 
ponential line that they may be represented:by. a ... 
single correlating cxirve. , This correlating line . ..:. 
is shown dotted in Fig. 8 and is given by the 
equation 

VF(B > 40) = 0.20 B‘2/^ : (16) 

This line is used to extrapolate the calculations 
to the Buoyancy number of 350 at the engine accel- 
eration of "0.02 g." At a density ratio of 10.0, 
which is close to the engine value, it Is seen in ^ 
Fig. 8 that the fuel volume fraction is decreased 
due to buoyancy effect by a factor of about:4.0. 

This reduction in fuel volume fradtibti &lcjuld,. 
be investigated and verified ih future flow 
lation experments.' The Buoyancy .number 6f :yari--, v 
ous engine designs will probably range from. ahpUt'dd 
B = 100 to 500, as shown in Fig. :8. . :By comparison^ 
present experiments operate at Buoyancy , numbers no,.c 
larger than about B = 30. Therefore h for Cacdh- :// 
rate assessment of the buoyancy effeef ,. tile : iahg^ 
of experimental Buoyancy numbers sho'dl-d be^ ex- 
tended into the range of engine values. 

Conclusions 


Effect of Acceleration on Fuel Volume Fraction 


Solutions of the flow equations were ob- 
tained for various fuel-to-propellant density ra- 
tios and for various vehicle accelerations. The 
vehicle acceleration was originsilly specified in 
terms of a Froude number, u|/gD, which is the con- 
ventional measui'e of fuel inertia to acceleration 
force. Values of Froude number were varied from 
10“^ to infinity. A simplified analysis indicated 
that density ratio should be included in a buoy- 
ancy parameter. Several dimensionless combina- 
tions are possible. After several different ways 
were used for plotting the calculated results, it 
was concluded that the best parameter for char- 
acterizing and generalizing the buoyancy effect has 
the form: 


A two-fluid steady Navier-Stokes equation so- 
lution has been programmed for the case of coaxial 
flow in a gas-core nuclear rocket cavity, includ- 
ing the buoyancy force due to density differences 
and vehicle acceleration. Solutions were obtained 
and fuel volume fractions were calculated for a 
propellant-to-fuel flow ratio of 50, a Reynolds 
number of 1000, in an L/D = 1 cavity. TTie fuel- 
to-propellant density ratio ranged from 1.0 to lo.o 
and the vehicle acceleration was varied from zcrit 
to a value for the engine. Calculated due., 
ume fractions are considerably below rsCues me. 
ttred in a coaxial flow experiment, ihis ir tr, ib- 
ably because the numerical soluticn convert/ s nip 
for Reynolds numbers below about iOOu ilicrt^u i ic 
experiment was at a Reynolds numbe t od iC". 

This report is concerned with the tr: lo • i'< hr 
calculated fuel volume fraction v:lu, , j : Llj i i i,m 
with the exact values. The Buoyancy 'innoc (iel .i'- 
mined in this report to best correiar- l.lu I'nol 


5 



volimie fraction with density ratio, fuel velocity, 
and vehicle acceleration has the form; 


(pp - Pp)gD 

PpUf 

Specific conclusions of this study are; 

1. At Buoyancy numbers below about 5, the fuel 
volume fraction depends strongly on density ratio, 
but is independent of vehicle acceleration. 

2. At Buoyancy numbers larger than 40, the 
fuel volume fraction is nearly independent of den- 
sity ratio, and decreases with the -2/3 power of 

Buoyancy nmber. 

3. For a typical engine flow with a density 
ratio of 10.0 and Buoyancy mmber of 350, the fuel 
volume fraction that would exist at zero g is de- 
creased by a factor of about 4.0 due to buoyancy 
effect. 

4. The range of current flow experiments 
should be extended to higher Buoyancy numbers to 
verify the calculated relationship between fuel 
vol:ume fraction, density ratio, fuel velocity, and 
Buoyancy number. 
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TABLE I. - VALUES OF THE GENERALIZED COEFFICIENTS IN EQUATION (l) 
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Figure 1. - Conceptual gas core nuclear rocket engine. 
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pp/pp= 1.0, 4.7 
mp/rrip >50 

L/D = iro ’ 

B = 0 - 30 
Re = 2x10^ 

VF <0.34 


(A) COAXIAL TESTS REF. (2). 



pp/pp ' 1. 0 10 

mp/mp = 50 

UD = 1.0 
B = 0 - 100 
Re = 1000 
VF< 0.072 


(B) ANALYTICALMODEL. 

Figure 2. - Comparison of coaxial flow con- 
ditions and analytical model. 



Figure 3. - Cavity model showing inlet profiles and buoyancy 
force. 
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(B) ppfpp = 4.0. 


Figure 4. - Streamline plots for different density ratios, 
mp/mp= 50, Re= 1000, B = 0. 



(B) pp/pp = 4.0. 


Figure 5. - Fuel mass fraction contours for different den- 
sity ratios, mp/mp= 50, Re= 1000, B = 0. 
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(A) B = 0. 



(B)B = 65. 

Figure 6. - Fuel mass fraction contours for different vehicle 
accelerations, mp/mp= 50, Re= 1000, pp/pp = 2.0. 



FUEL-TO-PROPELLAm- DENSITY RATIO, pp/pp 

Figure 7. - Variation of fuel volume fraction with fuel-to- 
propellant density ratio for zero acceleration. 
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Figure 8. - Fuel volume fraction as function of buoy- 
ancy number for various density ratios, mp/mp = 50, 
Re = 1000. 
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